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Abstract 

On the basis of Braaten and Segel's representation of the electro- 
magnetic dispersion relations in a QED plasma we check the numerical 
accuracy of several published analytic approximations to the plasma 
neutrino emission rates. As we find none of them satisfactory we de- 
rive a new analytic approximation which is accurate to within 4% 
where the plasma process dominates. The correct emission rates in 
the parameter regime relevant for the red giant branch in globular 
clusters are larger by about 10 — 20% than those of previous stellar 
evolution calculations. Therefore, the core mass of red giants at the 
He flash is larger by about 0.005A4© or 1% than previously thought. 
Our bounds on neutrino magnetic dipole moments remain virtually 
unchanged. 
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1 Introduction 



In stars with densities below nuclear, neutrinos are emitted by the plasma 
process 7 — > Vu, the photo process je~ — > e~Vu, the pair process e~e + — > z/z/, 
and by bremsstrahlung e~(Ze) — > (Ze)e~Vv. In Fig. 1 we show the regions 
of density and temperature where each of these processes dominates over the 
others. The plasma process is particularly interesting because it does not 
occur in vacuum and yet dominates the stellar emission rates in a large range 
of temperatures and densities. Notably white dwarfs and the degenerate cores 
of low-mass red giants fall in this parameter region. Evolutionary calculations 
for these stars can be tested with great statistical significance. For example, 
globular cluster observations allow one to determine the core mass at the 
helium flash to within about 0.012 A4 Q or 2.5% at a la statistical confidence 
level (Raffelt 1990, Raffelt and Weiss 1992). Therefore, it is important to 
identify possible systematic effects that would enter at this or at a larger 
level. If the standard neutrino loss rates are multiplied with a factor F v 
where F v = 1 refers to the standard case, the core mass at the helium flash 
varies approximately as 5M. C = 0.020 M. Q 5F U . This strong sensitivity to 
the neutrino luminosity has been used to constrain possible non-standard 
contributions such as the plasmon decay by virtue of neutrino dipole moments 
(Raffelt 1990, Raffelt and Weiss 1992). 

In practical stellar evolution calculations the neutrino loss rates are imple- 
mented in the form of simple analytic approximations. Widely used versions 
are those of Beaudet, Petrosian, and Salpeter (1967), Munakata, Kohyama, 
and Itoh (1985), Schinder et al. (1987), and Itoh et al. (1989). With regard 
to the plasma process, the results of Beaudet et al., Munakata et al., and 
Schinder et al. agree with each other to better than 1% if the same neutrino 
coupling strength to electrons is used, while Itoh et al. made an attempt to 
improve the accuracy of the rates at low temperatures. 

It turns out that all of these rates are relatively poor approximations 
for T ^ 10 s K which is relevant for low-mass stars. They were optimized 
for higher temperatures and correspondingly higher densities in the diagonal 
band of Fig. 1 where the plasma process dominates. At higher densities, 
however, the approximate photon dispersion relation that had been used in 
all of these works is a rather poor approximation (Braaten 1991) with the 
result that the above plasma emission rates are bad approximations every- 
where. In response, Itoh et al. (1992) have derived a new analytic formula 
which is claimed to fit the exact results to better than 5% in the region 
where the plasma process dominates. Independently, Blinnikov and Dunina- 
Barkovskaya (1992) have published rates which were optimized for low-mass 
stars. 



2 



Since then Braaten and Segel (1993) have devised an approach to the 
photon dispersion relation which, for given T and p, reduces the calcu- 
lation of the plasma emission rate essentially to the numerical evaluation 
of a few one-dimensional integrals. Their method for the first time offers 
a simple and practical approach to check the accuracy of all of the above 
emission rates as well as the approximations that were used to study non- 
standard neutrino emission (Raffelt 1990, Raffelt and Weiss 1992, Blinnikov 
and Dunina-Barkovskaya 1992, Castellani and Degl'Innocenti 1992). On the 
basis of this method we found none of the above approximations satisfactory 
and thus have derived a new analytic approximation to the standard rate 
for the plasma process which is accurate to better than 4%. For the non- 
standard emission rates we found that a simple scaling of the standard rates 
as in Raffelt and Weiss (1992) introduces only a small error (less than 5%) 
for low-mass stars so that bounds on neutrino dipole moments are mostly 
affected by the accuracy of the standard emission rates. 

We begin in Sect. 2 by adapting the results of Braaten and Segel (1993) to 
our calculation of standard and non-standard emission rates by the plasma 
process. Besides the standard neutrino interactions we include the possibility 
of the coupling of right-handed neutrinos to electrons (Fukugita and Yanagida 
1990), neutrino dipole moments, and neutrino electric "millicharges" . In 
Sect. 3 we perform a numerical test of all of the above analytic approxima- 
tions, and we derive and test a new version. We also check the accuracy 
of a simple scaling of the standard rates to obtain the non-standard ones. 
In Sect. 4 we discuss the accuracy of previous calculations of the core- mass 
at the helium flash as well as the accuracy of previous bounds on neutrino 
properties, and we derive bounds on the Fukugita and Yanagida model as 
well as on neutrino millicharges. 

2 Exact Emission Rates 

2.1 Electromagnetic Excitations in a Medium 

The behavior of electromagnetic excitations in a linear medium is best un- 
derstood on the basis of the wave equation for the vector potential A in 
Feynman gauge, [K 2 — II] A = j cxt , where j ext is an external electromagnetic 
current and II is the polarization tensor. Moreover, K 2 = uo 2 — k 2 where u 
is the frequency while k — |k| the modulus of the momentum coordinate. In 
an isotropic medium which also remains invariant under a parity transforma- 
tion, the polarization tensor can be expressed as (Nieves and Pal 1989a,b) 
II = 71tQt + t^lQl where Qt is a ^-dependent projector on the subspace 
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transverse to k while Ql projects on the longitudinal subspace, and in addi- 
tion Qt,lK = so that UK = as required by gauge invariance. Explicit 
expressions for Qt,l hi terms of K were given, for example, by Weldon (1982). 
He also showed that in the rest-frame of the medium ti^ = —U 00 K 2 /k 2 and 
vr T = - n L )/2. 

If II is expressed in this way, it is clear that the homogeneous wave equa- 
tion (K 2 — II) A = has non-trivial solutions only for 



UJ 
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k 2 — 7r s (u, k) = 0, s = T or L, (1) 



a relationship which implicitly defines the dispersion relations for transverse 
and longitudinal propagating modes ("photons" and "plasmons" or "trans- 
verse and longitudinal plasmons"). 

In order to calculate the decay rate into neutrinos of such modes we need 
to normalize properly the amplitude of the quantized excitations. To this 
end we consider a mode k, polarization s = T or L, with the corresponding 
frequency u, % k so that the dispersion relation Eq. ([]]) is obeyed. Then we 
expand it s (uj, k) = n s (u s ^, k) + w 27r s (a; S) fc, k){uj 2 — u 2 k ) so that the equation 
of motion for the amplitude A s ^ is that of a harmonic oscillator. Equivalently, 
we write the photon propagator near its poles in the form Z s /(oj 2 — ui 2 k ). 
Either way, one finds that the squared matrix element of a process with an 
external electromagnetic excitation of momentum k and polarization s carries 
a factor 

Z a =[l- <9 w 2 7r s (u; Sife , k)]' 1 . (2) 

In a fully ionized plasma the polarization tensor is simply given in terms of 
the forward scattering amplitudes of photons on the electrons and positrons. 
To lowest order in a = e 2 /4tt it is found to be 



fl^ = Sua 



d 3 p f(E p ) P-K (K^P U + K V P„) - K 2 P^P V - (P-K) 2 9flu 



(2vr) 3 E p (P-K) 2 -K a /A 

(3) 

where P = (E p , p) is an e or e + four-momentum. The sum of the phase- 
space occupation numbers for these particles is 

f( E ) = e (E-p)/T + I + e {E+rf/T + I ( 4 ) 

with \l the electron chemical potential and T the temperature. On the basis 
of Eq. (|3]) it is, in principle, straightforward to derive ir s and Z s for a given 
\i and T . 



4 



2.2 Analytic Representation for II and Z 

Braaten and Segel (1993) have shown that ignoring the K 4 /4 term in the 
denominator of Eq. (H) introduces an error which appears only on the a 2 
level so that to 0(a) it can be ignored. In fact, ignoring this term provides 
a better approximation to the 0(a) result than keeping it because ttt then 
remains real for all conditions as it must because the imaginary part from 
the 7t — > e + e~ process which otherwise appears at sufficiently large plasma 
frequencies is unphysical (Braaten 1991). 

Once the K A /4 term has been dropped, the angular part of the integral 
in Eq. (0) can be done analytically, leaving one with a one- dimensional inte- 
gral over electron momenta which can be done analytically in the classical, 
degenerate, and relativistic limit. In these cases Braaten and Segel foundP] 



TT T (U, k) 

7c L (u,k) 



1 + \G(vlk 2 /u- 
1 - G{vlk 2 /u 2 ) 



+ v 2 k 2 - 



(5) 



where t>* = ui/uip has the interpretation of a typical velocity of the electrons 
in the medium. The plasma frequency up and the frequency uj\ are defined 
by 



w P = 



Wi = 



4a 

7T 

4a 



7T JO 



dp (v 

i 

dp ( 



pf(E p ). 
r--v*)pf{E p ) 



(6) 



where v = p/E p is the electron or positron velocity. In the classical limit one 
has v* = (5T/m e ) 1//2 , in the degenerate limit v* = vp (velocity at the Fermi 
surface), and in the relativistic limit v* — 1. The function G is defined by 



G(x) 



3 

x 



2x 

y 



2v^ 



x 1 
log- 



6E 



x 



\ (2n + l)(2n + 3) ' 



(7) 



We have plotted G(x) in Fig. 2; note that G(0) = 0, G(l) = 1, and G'(l) = 
oo. 

Braaten and Segel then claim (and we have checked) that these results 
apply approximately for all conditions. The deviations between 7T s (uj, k) given 
by Eq. (|J) and by the proper integral over the e + e~ phase space are always 

1 Braaten and Segel used Coulomb rather than Feynman gauge so that we had to 
translate their longitudinal expression accordingly. We preferred to follow Weldon (1982) 
in the choice of gauge because the dispersion relations and plasmon decay rates then have 
the same form for transverse and longitudinal excitations. 
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so small that Eq. (|5|) can be considered exact to 0(a). As a higher precision 
would require an 0(a 2 ) QED calculation, nothing is gained by evaluating 
the full integrals. Therefore, the above results are as exact as an 0(a) result 
can be. 

In order to calculate the plasmon decay rates we also need the ampli- 
tude normalization factors. Braaten and Segel (1993) found the analytic 
representations 



^T,k 



J L,k 



^T,k( U T,k ~ V l k2 ) 



c4 )fe (3u;I - 2n T , k ) + + fc 2 )K >fe - t£*») ' 
2<,« fc -^ 2 ) 



[3^P - (w£, fc - ^fc 2 )]7TL,fc 



In these expressions 7r Sjk = ^(w^j;, fc) is the "on shell" value of 7r s for excita- 
tions with momentum k, i.e. u 2 sk — k 2 = n s>k . 

The decay 7 S — ► Vv is only possible if 7 S has a time-like four-momentum, 
^s,k > k. For longitudinal excitations this condition is only fulfilled for 
k < k mgx where 

n l/2 

(9) 



Up 



3 / 1 . l + <u* 
— ( — log 



2u. 



to the same level of approximation. 

In the Braaten and Segel representation the on-shell values ir S)k and thus 
the dispersion relations depend only on the medium parameters up and i>*. 
However, media are more readily characterized by their density and temper- 
ature. Therefore, in Fig. 3 we show contours of constant and 7 = Up/T in 
the region of p and T relevant for the plasma process in stars. In the shaded 
area the plasma process contributes more than 10% to the total neutrino 
luminosity. Evidently, it is important only for 0.3 ^ 7 ^ 30. Moreover, the 

contours are almost perfect vertical straight lines in this regime, i.e. the 
medium is degenerate. 

The polarization functions Tr s ,k/up and the amplitude factors Z S)k are 
functions only of t>* and of k/up. In Fig. 4 we show contours for these 
functions in the t^-fc-plane for transverse excitations. We have 1 < 7TT,k/wp < 
3/2 and Z?.k < 1- The deviation of Zx,k from unity is always small. In Fig. 5 
we show similar contours for the longitudinal case. We have n L)k /up < 1; 
the occurrence of the plasma process in addition requires < ~Kh,k- The 
contour ttl^ = is identical with the function fc max (f*) given in Eq. (|9|). The 
amplitude function Z^ k diverges when k — ► k max . Therefore, we show in 
Fig. 5 (lower panel) contours for Z* Lk = Z L ^ Ltk ju\^ k instead.^ 

2 Our Z* L k is what Braaten and Segel (1993) call Z L (k). 
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2.3 Plasmon Decay Rates 

In order to calculate the neutrino emission rates for the standard and several 
non-standard interaction models it remains to calculate the decay rates for 
the process 7 S — > Vv. In the standard model one finds (Adams, Ruderman 
and Woo 1963, Zaidi 1965, Dicus 1972) 



y2 

487T 2 « UJ S 



The effective vector coupling constant is 

C 2 v = jz C 2 v>i = (| + 2 sin 2 9w) 2 + 2(| - 2 sin 2 6 W ) 2 , (11) 
i=i 

where Cy^ is the effective neutral-current vector coupling constant of neutrino 
flavor i to the electrons. With a weak mixing angle of sin 2 = 0.2325±0.0008 
this is C v = (0.9312 ± 0.0031) + 2 (0.0012 ± 0.0001) = 0.9325 ± 0.0033 where 
the first term is from 7 S — > V e v e while the second term is from v^v^ and 
V T v T . Thus, even if these latter flavors were too heavy to be produced by 
plasmon decay the value of C v would change by less than its uncertainty. 
The contribution of the axial neutral current is always negligible. 

As a first non-standard coupling we consider a model by Fukugita and 
Yanagida (1990) which was devised to obtain a large neutrino magnetic dipole 
moment and a strong effective coupling of right-handed u e 's to electrons. The 
relevant part of the Lagrangian is 

C VReL = g<N>u R A L + he. (12) 

where is a scalar field of mass M and g is a dimensionless coupling con- 
stant. For low energies this interaction produces an effective neutral-current 
coupling of r.h. neutrinos to electrons. The corresponding effective vector 
coupling constant is 

C ^ R = 32M*G 2 • (13) 

In this model the rate for 7 S — > v e ^v e ^ is found by inserting C VR instead of 
C v into Eq. ©. 

Next, we consider direct couplings of neutrinos with the electromagnetic 
field. The least exotic possibility is that of neutrino dipole moments with an 
effective Lagrangian 

1 3 _ __ 

^ = ^E (m&iP^j + ^^^vi^j) , (14) 
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where /i^- and are matrices of magnetic and electric dipole and transition 
moments, F is the electromagnetic field tensor, and the sum is over neutrino 
flavors. We define an effective dipole moment by 

A= E(>d 2 + M 2 ) ( 15 ) 

with the restriction that the sum should only run over those flavors which 
are light enough to be produced by the plasma process: mi + rrij % wp. Then 
we find for the plasmon decay rate 

r sfc = ^^k (16) 

This result is in agreement with Sutherland et al. (1976) except for their 

Finally, we assume that neutrino flavor i carries a "millicharge" q^e. In 
this case we find 

T-l ll a Z S,k^S,k / 17 N 

i-s,k~— , {l<) 

<J u s ,k 

where a = e 2 /4ii is the (electron) fine structure constant and ql = X)f =1 q 2 . 



2.4 Neutrino Emission Rates 

If transverse and longitudinal electromagnetic excitations of momentum k 
can decay according to 7 S — > Vv with a rate Y s j~ then the energy-loss rate 
per unit volume of a medium at temperature T is Q = Qt + Ql where 



Q 



T 



dk k 2 ry fe uj T ,k 



2n 2 e^/ T - 1 ' 



r k ^dkk 2 T L:k u L:k 
Ql ~ Jo ~2^~ ~ i • ( 18 ) 

Inserting the results of the previous section into this equation we find for the 
various interaction models 

G 2 

Qv = C v T ujpQ 3 , 

Q q = |^TVQi- (19) 
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The dimensionless emission rates are Q n = (Q^ + Qn) where 



Ql = 2/ 

JO 



00 dkk 2 _ ( ir Tik \ n 1 



T 3 ^ V u 2 ) e^/T -1 ' 

In the Braaten and Segel approximation the Q n are only functions of v* and 
7 = up/T. 

In Fig. 6 we show contours in the iv7-plane of Q3/Q3, i.e. the ratio 
between longitudinal and transverse emissivity for the standard-model in- 
teractions. Corresponding contours in the p-T-plane are shown as dashed 
lines in Fig. 1. Evidently, the longitudinal process is of importance only in a 
relatively narrow region near 7 = 10. 

In a practical calculation of anomalous neutrino losses one will scale the 
standard luminosity appropriately. The relevant ratios are 

Q, _ t42na Q 2 _ Q /l0keV\ 2 Q2 ^ ^ 

Qv CyGpujp Qs \ ujp J Q 3 ' 

where /ii 2 = fi u /l0~ 12 (e/2m e ) and 

where 514 = q u /10~ w . In Fig. 7 we show contours of Q1/Q3 and Q2jQz- 



3 Numerical Emission Rates 

Even though the methods described in Sect. 2 allow for a relatively simple 
calculation of the plasma neutrino emission rates, one still needs to evaluate 
several numerical integrals in order to determine the energy loss rate for given 
values of density and temperature so that this procedure is not suitable to 
be coupled directly with a stellar evolution code. However, we can use these 
results to test the accuracy of widely used analytic approximation formulae. 

It turns out that for the plasma process the analytic approximations of 
Beaudet, Petrosian, and Salpeter (1967), Munakata, Kohyama, and Itoh 
(1985), and Schinder et al. (1987) all agree with each other to an astonish- 
ing accuracy if the same value for C v is used. In Fig. 8a we compare the 
numerical rates of Beaudet et al. (1967) with the exact results obtained by 
the Braaten and Segel method. We show contours for the relative deviation 
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in percent, Q analytic /Qcxlct — 1- We stress that what is plotted is the error 
of the total emission rate; the plasma process alone may be more inaccurate 
than shown in regions where it does not dominate. In this and the following 
figures we have used the rates of Itoh et al. (1989) for the other emission pro- 
cesses. Thus, Qanaiytic an d QcxLt differ only in the treatment of the plasma 
process. From Fig. 8a it is evident that the numerical rates are rather poor 
approximations almost everywhere. Nevertheless, it is these rates that have 
been used in virtually all stellar evolution calculations. 

Recently Itoh et al. (1992) have published numerical rates for the plasma 
process which are claimed to be more accurate than 5% wherever the plasma 
process dominates. In Fig. 8b we put this claim to a test and find that there 
remain substantial regions with much larger errors. We have checked that 
at least some of these deviations also occur between their tabulated emission 
rates and their fitting formula. 

Blinnikov and Dunina-Barkovskaya (1992) have also derived a new ana- 
lytic approximation which is optimized in the region of small temperatures 
and densities relevant for low-mass stars. We compare their rates with the 
exact ones in Fig. 8c. Indeed, these rates are rather good fits for T ^ 10 s 
and p ^ 10 6 , but for higher T or p they are unrelated to the exact results. 
These approximation formulae involve multiplicative factors which depend 
on T alone and which approach 1 for T — > 0. If we set these factors to 1 
for all T the errors of the resulting simplified emission rates are shown in 
Fig. 8d. The fit is not worse in the low-T and low-p region, and much better 
otherwise! 

As we find none of the published fitting formulae satisfactory we have 
derived a new one. To this end we have started with the T = version 
of Blinnikov and Dunina-Barkovskaya (1992) and then "flattened" the errors 
with an extra correction factor f xy . As a result we have come up with Qpias = 
CyQapprox (in ergcur 3 s _1 ) where 

Qapprox = 3.00 x 10 21 A 9 7 6 e" 7 (f T + f L ) fxy , (23) 

where A = T/m e and 7 = ujq/T where u is the zero-temperature plasma 
frequency, Uq = 4naN e / E F . Numerically, 

A = 1.686xlO _10 T 

1.1095x10^/^ 

7 T 2 [l + (1.019X10- 6 p//v) 2/3 ] 1/2 ' 1 ' 

with T in K, p in g/cm 3 , and p e the number of baryons per electron. More- 
over, we have 

f T = 2.4 + 0.6 7 1/2 + 0.5l7+ 1.25 7 3/2 , 



10 



h = 



8.6 7 2 + 1.357 7/2 
225 - 177 + 7 2 



(25) 



The coefficients here are slightly different from those used by Blinnikov and 
Dunina-Barkovskaya (1992). Finally we define 



x 



1 

6 
1 

V = 6 



+ 17.5 + log 10 (2p//i e )-31og 10 (T) 
-24.5 + log 10 (2p//i e ) + 31og 10 (T) 



(26) 



where x is a coordinate transverse to the diagonal band in Fig. 1 where the 
plasma process is important, and y is along this band. If|x|>0.7ory<0 
we use f xy = 1 and otherwise 



fxy = 1-05 + 



0.39- 1.25 x- 0.35 sin(4.5a;) - 0.3 exp{-(4.5 x + 0.9) 2 } 

2 



X 



x exp 



min(0, y — 1.6 + 1.25 a;) 
0.57 - 0.25 a; 



(27) 



We show the errors of our fitting formula in Fig. 8e; it is better than 5% 
everywhere and better than 4% almost everywhere. 



4 Discussion and Summary 

In low-mass stars helium ignites in the cores of red giant stars under de- 
generate conditions. Helium burning depends extremely sensitively on tem- 
perature and density so that even relatively minor changes in the neutrino 
cooling rates of the core change the ignition point and thus the core mass 
Aitip at the tip of the red giant branch, which in turn affects the luminosity 
during the subsequent horizontal branch evolution. All previous calculations 
of the core mass at helium flash seem to have used the Beaudet et al. (1967), 
the Munakata et al. (1985), or the Schinder et al. (1987) rates, all of which 
are equivalent with regard to the plasma process and thus underestimate the 
emission rate as shown in Fig. 8a. 

In order to illustrate the relevant range of parameters we show in Fig. 9 
the evolution of the central density and temperature of a 0.80 A4 Q star where 
the tail ends of the arrows mark the conditions when the surface brightness 
is at the indicated magnitudes. In Fig. 10 we show the red giant part of 
this track overlaid with the errors of the Munakata et al. (1985) rates that 
were used in our previous calculations (Raffelt and Weiss 1992). The error 
contours are virtually the same as those shown in Fig. 8a. Near the helium 
flash the average neutrino emission rate was underestimated by around 15%. 
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Raffelt and Weiss (1992) found that the core mass at the helium flash 
varies approximately as SA4 C = 0.0205F U if the standard neutrino loss rates 
are multiplied with a factor F u , a result which agrees with the previous 
calculations of Sweigart and Gross (1978). Therefore, the core mass at the 
helium flash will be increased by about 0.004 .M©. In order to confirm this 
estimate we have re-calculated run 11 of Raffelt and Weiss (1992) with the 
analytic emission rates derived in this paper. For M. = 0.80.M©, Z = 10~ 4 , 
and Yq = 0.22 we previously found Ai t i P = O.497A^ for the core mass 
at helium ignition (at the tip of the red giant branch) while we now find 
Map — 0.503.M©. Therefore, we find an increase of <LM t i P = O.OO67W , 
in reasonable agreement with our simple estimate. Thus, using the correct 
neutrino emission rates changes the core mass at the helium flash by a small 
but noticable amount. 

In order to constrain neutrino dipole moments Raffelt and Weiss (1992) 
as well as Blinnikov and Dunina-Barkovskaya (1992) and Castellani and 
Degl'Innocenti (1992) have included the increased plasma losses in evolu- 
tionary calculations by scaling the standard rates with a certain factor. 
Instead of the exact ratio given by Eq. ([H]) they used Q2/Q3 = 1 and 
Up = Uq = AiraN e / Ep which is the zero-temperature value for the plasma 
frequency. In Fig. 11 we show the evolutionary track of Fig. 9 in the plane 
of 7 = up/T and f*. Evidently on the RGB the interior of the core has an 
almost constant value 7 pa 3. A comparison with the upper panel of Fig. 7 
reveals that by using Q2/Q3 = 1 Raffelt and Weiss (1992) have underesti- 
mated the dipole-induced emission rate by about 5%. An additional small 
error occurs by using u Q instead of up. In Fig. 12 (upper panel) we show 
the error of the dipole-induced emission rate if it had been scaled with the 
correct standard rate. These errors are so small that this scaling procedure 
remains well justified. 

In Fig. 12 (lower panel) we show the compound error of the dipole-induced 
emission rate when coupled with the Munakata et al. (1985) emission rates. 
Near helium ignition the emission rate was underestimated by 15 — 20%. 
Because the core-mass and brightness at the helium flash vary approximately 
linearly with \x v in the range of interest, for a given value of \i v these quantities 
are changed by about 10% more than given in Raffelt and Weiss (1992). This 
is a negligible change with regard to bounds on \l v . 

The bounds on neutrino dipole moments discussed in Raffelt (1990) and 
Raffelt and Weiss (1992) crudely amount to the constraint that the total 
neutrino luminosity must not exceed about twice its standard value in order 
to maintain the beautiful agreement between the observed and calculated 
properties of globular cluster stars. We have already discussed that on the 
RGB we may set Q2/Q3 = 1 in Eq. (|2"T| ) without introducing a large error, 
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and similarly we may set Q1/Q3 = 1 in Eq. fl22|), although the error here 
is somewhat larger. In any case, these approximations lead to an underesti- 
mation of the non-standard emission rates and thus to conservative bounds. 
Moreover, near the helium flash the density is about 10 6 g/cm 3 (see Fig. 10) 
so that Up = 17.8 keV in the center of the star. Therefore, in the center of 
the star we have for the total energy-loss rate from Eqs. (|2TD and (f22|) 



Away from the center the density and thus uj? is lower so that the coefficients 
of n\ 2 an d of ^14 would be larger if averaged properly over the core. Hence 
this expression, again, is a conservative estimate of the non-standard neutrino 
losses. With Qtot/Qv < 2 as a formal criterion we find the bounds 



In the model of Fukugita and Yanagida (1990) one has both a dipole moment 
and r.h. interactions for the z/ e 's so that in that model the constraints on the 
dipole moment are more restrictive. The bound on [i v alone is virtually the 
same as that from the more detailed analysis of Raffelt (1990) and Raffelt 
and Weiss (1992). The bound on q v is slightly more restrictive than that 
found by Bernstein, Ruderman, and Feinberg (1963). 

In summary, we have discussed in detail the neutrino emission rates from 
the plasma process due to standard and non-standard interactions. By means 
of Braaten and Segel's (1993) representation of the QED dispersion relations 
we have tested the accuracy of widely used analytic approximation formulae, 
none of which are found satisfactory. For the first time we have derived an 
approximation formula which is accurate on the 4% level wherever the plasma 
process dominates. The correct emission rate leads to a slightly increased 
core mass at the helium flash in low-mass stars, and to a slightly increased 
sensitivity to non-standard neutrino losses. While these changes are noticable 
they are so small that previous bounds on neutrino dipole moments remain 
virtually unchanged. 



Qtot/Qv = 1 + 1-07 C V>R + 0.100 /4 + 0.066 q 2 u . 



(28) 



C v ,r < 0.9, 
\i v < 3xl0~ 12 (e/2m e ) 
q u < 4xl0~ 14 . 



(29) 
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Figure Captions 



Figure 1 

Regions of density and temperature where the different neutrino emission 
processes contribute more than 90% of the total. \x e is the mean number 
of baryons per electron. The bremsstrahlung contribution is for helium. 
The dashed lines are contours for the indicated values of Ql/Qt, i.e. the 
contribution of the longitudinal relative to the transverse plasma process. 

Figure 2 

Function G(x) as defined in Eq. (0). 
Figure 3 

Contours of 7 = Wp/T and t>* as defined in Eq. (H). 
Figure 4 

Contours for n Tjk = uj\ k — k 2 in units of Up and for Z Tfk . 
Figure 5 

Contours for n L ^ k = u 2 Lk — k 2 in units of co>p and for Z* Lk = Z L)k / K L ^ k /uj\ k . 
The n L)k = contour corresponds to fcmax^*) of Eq. (||). In the cross- 
hatched region longitudinal plasmons have a space-like four momentum and 
thus cannot decay. 

Figure 6 

Contours for Q3/Q3 as defined in Eq. (p0|). This ratio is identical to the ratio 
of the longitudinal and transverse emission rates in the standard model. 

Figure 7 

Contours for Q2/Q3 and Q1/Q3 as defined in Eq. (|20|) . 
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Figure 8 

(a)-(e) Comparison of the analytic plasma emission rates of the indicated 
authors with the exact rates. Except for (e) the contours are at multiples of 
±10%. The errors in (a) of the Beaudet et al. (1967) rates are the same for 
Munakata et al. (1985) and Schinder et al. (1987). 

Figure 9 

Evolution of the central density and temperature of a 0.80 Ai Q star up the 
red giant branch (RGB), and then from the horizontal branch (HB) up the 
asymptotic giant branch (AGB). The rear-ends of the arrows are at the in- 
dicated values for the surface brightness. 

Figure 10 

Error in % of the standard neutrino emission rates used in the evolutionary 
calculations of Raffelt and Weiss (1992). The evolutionary track is that of 
Fig. 9. 

Figure 11 

Evolutionary track of Fig. 9 in the 7-v*-plane. 
Figure 12 

Error in % of the dipole moment emission rates used in the calculations of 
Raffelt and Weiss (1992). Upper panel: Error relative to the standard rates. 
Lower panel: Compound error if combined with the analytic standard rates 
of Munakata et al. (1985). 



16 



